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Abstract 

The Bohr-Sommerfeld rule for a spin system is obtained, including the first quantum corrections. 
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model, and found to agree very well with exact results. 
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The question in this Letter was investigated twentysome years ago by R. Shankar, in a 
paper in this journal with a nearly identical title Many physical systems involve large- 
magnitude spin or pseudospin degrees of freedom. Examples include spinning molecules 
3, 3], superdeformed rotating nuclei P], magnetic molecular solids and the Lipkin- 
Meshkov-Glick (LMG) model of certain collective excitations in nuclei p]. In all these cases, 
because the spin magnitude, j, is large, a semiclassical approach is natural. In some cases, 
such an approach is effective even for j as small as 1/2, e.g., in studies of magnetic order 

It may therefore surprise some readers that so far we have lacked even the Bohr- 
Sommerfeld (BS) quantization rule for a single spin. All previous authors cited have used 
leuristic or ad hoc BS rules that vary in their treatment of the first quantum corrections 
lo| . In this paper, we use the semiclassical propagator for spin, which has only recently 



been fully understood to find the BS rule systematically. We shall first state the rule, 
then apply it to a few model problems, and finally prove it. 

The general problem is as follows. We are given a quantum Hamiltonian 

n = niJ:^,j°^,r/), (i) 

where {i = x, y, z) are components of the dimensionless smn operator, J°'', with the 
commutator [J°^ , J°''] = ieijkJ'^^, and 7"^ is a polynomial in J°'' jl^]- Since J°p- J°p = j(j + l) 
is a constant of motion, the spin magnitude j must also be given. We want the eigenvalues 
of H when j is large, and the problem is, intuitively, quasiclassical. 

The first step is to understand the classical limit. We expect Ti. to correspond to some 
classical Hamiltonian Hc{J) of a c-number angular momentum J, with dynamics defined by 
the Poisson brackets {Ji, Jj} = €ijkJk- Since J ■ J is conserved, we can regard the motion 
as taking place on a sphere of radius j = |J|. All orbits on this sphere are closed and non 
intersecting, so it may regarded as the phase space of the system. With polar coordinates 
6 and (p, and the usual identification of and (p as canonically conjugate momentum and 
coordinate, we can expect the BS rule to quantize the action integral / cos 6d(j). Robbins et al. 

set this integral to (2n+l)7r/ + Barter and Patterson Q set it to 2n7r/[j(j+l)]^/^, and 
Shankar [l| equates jA^^ [see paragraph below Eq. (j27j) ] to 2n7r. Neither rule is formulated 
in a context that extends to half-integer j. The rule we find is 

1 fTiE) 

(j + l)At{E) + -J^ J VjH^dt = {2n + l)7r. (2) 



Here, n is an integer ranging from to 2j, and 



At{E)=i (l-cos^)rf0 (mod47r) (3) 

is the area enclosed by an orbit of energy E jl^ . For future use, let us denote the two terms 
in Eq. (j2I) by Ji and I2, and also define 



(4) 



The term I2 is also evaluated along the orbit He = E, which is taken to have a period T{E). 
This term, and the extra \ in the j + \ factor in Ii represent the first quantum corrections. 

To make Eq. Q precise, we must also state the rule for associating with Ti. We do 
this by writing 

(5) 



£=0 m=-j 



where are the spherical harmonic tensor operators, and the coefficients Qm are uniquely 
determined since Tr(3^|^3^£/m/) oc 5«'5mm'- Then, 



(6) 



with the same Qm's, and the Y^m's are solid harmonics. Along with the J^^m's, they may be 
defined by the generating function [l^ 



(a-J)^ 
(a - J°P)^ 



Atx 



m=-t 



(7) 



where, a^m = ^!/[(^ + m)\{i - m)!]^/^ and 

. A 



1 



-(x + zy) + — (x-zyj 



(8) 



is a complex vector with a ■ a = 0. Note that given Ti, is completely determined. To 
get He from H, we also need to give the value of j. Hence, the H ^ H^ correspondence is 
one-to-one. Note also that the trinomial expansion of (a- J°p)^ is analogous to the expansion 
of {aq + bp)"- that gives Weyl ordering. 

Let us now apply Eq. to a few models. Our first example is = ujJ"^, and Hc = ojJz- 
The orbits are given by cos^ = E/ujj = const, and (f) = uj. Hence, Ji = 2j7r(l — E/ju), and 
I2 = nE/ju!, so the quantization condition is En = uj{j — 11)11, with n = 0, 1, 2, . . ., 2j. In 
this simple case, the rule is exact. 



Our next example is somewhat less trivial, 7i = z/(J°p)^, He = z/J|. Again cos^ = const, 
E = vj"^ cos^ 9, and (f) = 2uj cos 6*. The rule again yields j cos6 = {j — n), i.e.. En = v{j — ny'. 
This is also exact, but more importantly, it is in accord with Kramers' theorem. It is not 
hard to see that this is true for any time-reversal invariant 7i, for which Hc{—J) = Hc{J). 

As our third example, we consider the LMG model, just as done by Shankar. The 
Hamiltonian now is 

and He is obtained by simply deleting the 'op' suffixes. On the orbit with energy E, 



^ 1 ± [1 + rcos20(rcos20- 2E/j)li/2 

cost^ = ■ — . (10) 

r cos 2(p 

Now, He —He under a 180° rotation about x ± y, i.e., O^n — 9,(j)—>-(j) + 7i/2. This 
symmetry forces states to occur in pairs, E and —E, plus a nondegenerate state at E' = 
for integer j. For orbits related by this symmetry, we have /i 4:71 j — Ji, and I2 —>■ so 
the BS spectrum is also symmetric about E = 0. For the orbit at = 0, Ji = (2j + l)7r, 
and /2 = 0, so = is an allowed energy for integer j only. From now on, we consider only 
E>0. 

The energy landscape for this model looks like this. If r < 1, i^c has a maximum at 
6 = {He = j) and a minimum a.t 6 = n. If r > 1, these points turn into saddle points, and 
He develops two degenerate maxima [with He = |j(r + r^^)] along the = and = vr 
meridians. Orbits at slightly lower energies circle these maxima, and are separated by a 
separatrix a.t E = j passing through the north pole. (See Fig. 1 in Ref. Q.) Hence for r > 1, 
we expect to have pairs of levels split by tunneling for E above or around j, and single levels 
below. In this paper, due to limited space, we will not incorportate tunneling effects into our 
BS calculations. In principle this may be done by allowing the orbits to become complex y | . 
We have discussed how to find ground pair splittings from the propagator elsewhere |l5|. 
For the orbits with E > j, both signs in Eq. (ITUI) are valid, but for E < j (and any r) only 
the minus sign is meaningful. 

The evaluation of the two action terms in Eq. Q requires simple one-dimensional numer- 
ical integration. The resulting BS spectrum is compared with the exact one (from numerical 
diagonalization of 7-^) in Table B The values of j and r are the same as those used by 
Shankar, and it is evident that our BS analysis improves on his. Indeed, it is rather good 
almost uniformly, but a few aspects call for comment. For r < 1, the highest energy exceeds 
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°By extrapolation 

'By mapping to p variables 



j, and can be found by extrapolating the BS action. This amounts to allowing for complex 
orbits. For r = 1, the energy is not quadratic in deviations about the maximum at = 0. In 
terms of the stationary phase approximation on which Eq. Q is premised (see below), this 
corresponds to an exceptional case where we must include fluctuations higher than second 
order (Gaussian) about the stationary phase point. The two highest energies in Table|l]were 
obtained via the mapping J°p ^ g, J°p ^ jp, with = and textbook BS quantization 
for a particle in one dimension. The same limitation on Eq. is present for r > 1 close 



5 



to the tunneling barrier, i.e., E ~ j. Now, one may have a pair of tunnel split levels such 
that one level is below the barrier and the other above. Literal use of Eq. (0) then yields 
two degenerate levels just below the barrier, with large errors. We have chosen to find the 
unbound partner by extrapolating the action ioi E < j. A proper semiclassical approach 
would require including complex orbits, and non-Gaussian fluctuations. These aspects are 
not unique to spin, and also occur with one-dimensional double well potentials. 

In the rest of the paper, we show how we derive the BS rule (j21). Let |n) be a spin 
coherent state with maximal spin projection along the direction n with polar coordinates 
(6', 0), i.e., J°P ■ n|n) = j|n). In terms of stereographic coordinates 



(11) 



the state may also be written as \z) = exp(zJ°'')|z). Note that \z) is not normalized, and 
{z\z) = {1 + zzYK We shall write |n), \z), or 16*, 0) interchangeably as needed. Secondly, we 
shall need to discuss points on the complex unit sphere, for which 9 and are not real, or 
equivalently, z and z are not true complex conjugates. Such points are specified by the pair 
{z,z). ^ 

Our starting point is the semiclassical approximation to the propagator K = 
{zf\e-'^^\z,y. 

1/2 / „■ ,T ■ 



Ksc{zf,Zi;T) = \ — 



2j \dzfdzi 



exp S + i^ Adt 



(12) 



with. 



S = j\nN + 



zz — zz 



1 + zz 



iHsciz,z) 



dt. 



A 
N 



+ z ^ z, 



d {1 + zzy dH, 
dz Aj dz 
{l + Zfz{T)){l + mz, 



(13) 

(14) 
(15) 



Here, S is the action along the least action trajectory from the point (^(0), Zi) to {zr, z(T)); 

I II 

as noted in Ref. 11, to find such a trajectory, i.e., a solution of the Euler-Lagrange equations. 



dz 



(16) 



we must allow ^(0) ^ zi and zj ^ z{T). The integral in S must be evaluated along this 
trajectory. So must / Adt, the Solari-Kochetov correction, which is 0(l/j) relative to S, and 



key to getting a propagator that is self-consistently replicating under composition. Further, 



H,,{-z,z) = {z\n\z), (17) 

which we shall call the semiclassical Hamiltonian. Finally, we must sum the the right hand 
side in Eq. (jl2j) over all solutions of the equations of motion if there is more than one. 

Our goal is to find the Green function, and then the energy spectrum by looking for its 
poles. In step 1, we evaluate the Laplace transform 

Fsc(%, z,- E) = r K,,{zf, zf, T)e'''^dT (18) 

•J 

by the stationary phase method. This naturally leads to the action at fixed energy W{E) = 
S{T) + iET, where E and T are related by 

dw(E) asm . . 

Then, using known methods Q], we can show that 



(20) 

T=T{E) 



Here we have included a Maslov-like phase e*"^, which arises because the mapping from T to 
E is many to one. This is best seen by considering the case zj = Zi, which we shall shortly 
encounter when we take the trace. Now, for the same energy E, there is more than one 
trajectory corresponding to multiple traverses of the same fundamental orbit. So, T is an 
integer multiple of the basic period Tq, the different branches of W differ by additive integer 
multiples of iETq, and the phase e^'^ is (—1)" for n traverses. 

In step 2, we again use stationary phase to perform the trace that gives the semiclassical 
Green function: 

^ 2? + 1 /■ dzdz „ , ^, ,„ . 

G^iE)^^l^^^-^F.A^.,.,E). (21) 

The stationary phase condition yields z{T) = z, and z{0) = z, so now we are only considering 
closed classical orbits with energy E, with momenta that match smoothly at the end points. 
Next, we must integrate over small fluctuations t] and f] in z and z. For the fundamental 
orbit at energy E, we can write 

= -JIJ eMZsc) exp(-lg) (22) 



using known classical mechanical identities. Here, 
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The quadratic form Q has a zero eigenvalue, corresponding to a fluctuation (f/, rj) that 
moves us along the orbit. If we parametrize this direction by the time t, and the orthogonal 
direction by a variable r, then. 



(25) 



and Q = 4r^|i;p. The integration over all points along the orbit turns into an integral over 
t from to T, and with the jacobian d{'r],f])/d{t,r) = — 2|i;p, the r integral may also be 
done. The result is that for one traverse of the orbit we get — iTexp(Zsc). For n traverses, 
we multiply by (—1)" and let Zg^ nZsc- Summing over all n, we get 

' l + expi-Z^E))- 
This has poles (each with unit residue) whenever 



(26) 



Zsc(^) = {2n + l)i7r. 



(27) 



Equation ()27|1 is our quantization condition. The first term in Eq. is —ij times the 
area A^^. Unlike A^, A'^^ is computed on the orbit H^c = E, not He = E, since it is Hsc 
that appears in Eq. (fTB|) . In the second term, we note that 



1 + zz) 



dzdz 



(28) 



where Vq is the angular part of the Laplacian. 

A true BS rule must be in terms of the classical Hamiltonian, He-, not Hsc- To cast 
Eq. ^TJ\ into the form Q, we must relate He and H^c- Using Eq. (|7j), we can show that 

- 1 



(n|3^,^(J°P)|n) 



with J = jfi |17|. Since, V^F£„(J) 
follows that i^sc = He + Hi, where 



1 - 



4j 



-£(£ + l)r,„(J), and J ■ Vjr,„(J) 



Hi = j,VlHe+^3-VjHe. 



(29) 

£y,^(j), it 

(30) 



We have written He instead of i^sc on the right since the corrections are of 0(j~^). 

The difference between ^+ and is now easy to find. Let us denote the solutions of 
Eq. (fTH|) (^sc5-2sc) and those with instead of ifgc by (^c^c)- If we regard and Zc as 
functions of E and then 

H,{z,,{E,z),z) = E-H,, (31) 
H,izciE,z),z) = E. (32) 

Writing Zgc = Zc + Zi, z^c = Zc + Zi, and using the equations of motion, we find 




In this way we find 

jAt=jAt+ r H,dt. (34) 

JO 

Substituting this in Eq. (j22I), and making use of Eqs. (PHj) and (jSHI), we see that Eqs. (|77j) 
and are equivalent. 

We mention in closing that we can also use coherent states to find a BS rule for one 
dimensional potential problems. The result is akin to Eq. (j2ZI)- To show its equivalence to 
the textbook rule, we use Weyl ordering and the Wigner-Moyal formalism to relate 7i, He, 
and Hsc- In the spin case, the problems of ordering, and of relating H^. to H^c-, are neatly 
solved by the generating formula (|7|). 
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